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We report a clear anomaly in the rf response of a tunnel-diode resonator measured in borocar- 
bide crystals in dc fields along the c— axis. We associate the anomaly in the dynamic magnetic 
susceptibility, x, with a structural transition in the vortex lattice, x sensitive to the Campbell 
penetration depth Ac of the rf perturbation into the mixed state. Ac depends on the vortex lattice 
elastic moduli, which in turn depend on the vortex lattice structure. The high-field transition line 
is clearly mapped and shown to behave in agreement with theory. 

PACS numbers: PACS: 74.25. Nf, 74.25. Qt, 74.25.Sv, 74.70.Dd 



In the classical picture of type-II superconductors, vor- 
tices should form a hexagonal lattice in isotropic or cubic 
materials; the same should happen in tetragonal materi- 
als for applied magnetic fields along the crystallographic 
c axis. It has been known for a long time that exper- 
imentally this is rarely the case [H, The reason for 
such behavior was, to a large extent, clarified after the 
nonlocal effects were incorporated into the London de- 
scription of intervortex interaction [1, 0, Q ■ This effort 
resulted in the realization that the vortex lattice (VL) 
structure depends on the Fermi surface, the order pa- 
rameter symmetry, as well as on the applied field, tem- 
perature, scattering strength, and even on thermal fluc- 
tuations of vortex positions ^0] . Currently, vortex lattices 
structural phase transitions are best studied in the rare 
earth-nickel-borocarbides 0, H, i, E, El, El, H [11 El , 
but in fact present in m any superconductors including 
cubic Nb or VaSi [il,E3,Ei. 

Hence, in addition to boundaries of different dynamic 
behavior (e.g., hquid, glass, or the "Abrikosov crystal"), 
the phase diagram of the vortex state contains lines sep- 
arating different periodic structures within the domain 
of the Abrikosov crystal. Information about these struc- 
tures is consequential because the VL elastic properties 
are structure dependent, whereas the dynamic behavior 
is affected by the lattice elastic moduli |19j . 

Vortex lattices have been probed by Bitter deco- 
ration [ISl, Small Angle Neutron Scattering (SANS) 
i EM Ho, E3| , and by Scanning Tunneling Microscopy 
111 . (Hereafter we are talking about tetragonal crys- 
tals in fields along the c-axis). Notwithstanding consider- 
able progress made in locating the transition lines in low 
fields (an example is given in Fig. [3]), the complete map- 
ping of transitions had proven elusive mostly due to poor 
high-field sensitivity of these techniques. Calculations 
based on the Ginzburg-Landau approach [ll[ suggested 
that the transition line for transformation of rhombic (tri- 
angular) to square (R-S) structures may end up at the 
upper critical field Hc2{T) . However the SANS data on 
LuNi2B2C have shown the transition line bending away 
of Hc2{T) [l5|. Gurevich and Kogan suggested that 



softening of the lattice near Hc2{T) causes the vortex 
fluctuations to increase thus smearing a weak square- 
symmetric intervortex interaction and making it isotropic 
(see also [29|). The model predicts that in small fields the 
R-S transition line Hu (T) slowly increases with tempera- 
ture but then sharply turns upward to avoid Hc2{T) and 
become a double- valued function. 

The dynamic magnetic susceptibility % of a supercon- 
ductor in the mixed state can be measured by the tunnel- 
diode resonator technique described elsewhere |2l|, i23 |. 
In short, the frequency shift. A/, due to a small change 
in the resonator induction caused by the finite penetra- 
tion of the ac field into the sample in the mixed state 
is recorded. This shift is proportional to the change in 
the total rf penetration depth, S oc A/ as long as S is 
smaller than the sample size. In the limit of small ampli- 
tude ac excitation fields (0.02 Oe in our case), the non- 
dissipative part Re x is determined by the ac field pene- 
tration depth S = (A^ -I- A^)^/^, where A is the London 
penetration depth and the Campbell length Ac charac- 
terizes the pinning and the VL elastic response 23|, [24 1 . 
For / ^ 10 MHz, the skin-depth due to the normal exci- 
tations vastly exceeds the size of our samples and has no 
effect on 5. In this work we utilize the sensitivity of this 
method to small changes in Ac that accompany the VL 
transitions to map the high- field part of the i/g (T) line 
unattainable by other techniques. 

LuNi2B2C (Lul221) were grown out of Ni2B flux 
[13, [H- Powder X-ray diffraction spectra indicate 
that there were no detectable second phases present. 
YNi2B2C was grown using vertical zone melting method 
with a commercial 4-mirrors image furnace (model FZ- 
T-4000-H-VLVPM-PC, Crystal Systems Corp., Japan). 
We have studied both as grown and annealed in ac- 
tive vacuum (950°C, 48 h) samples of a typical size 
1 X 1 X 0.1 mm'^. Residual resistance ratios (RRR) of 
Lu-1221 and Y-1221 were about 25, the critical temper- 
atures Tc « 16.2 and 15.2 K and the upper critical fields 
Hc2 ~ 9.2 and 8.8 T, respectively. Hc2 was determined 
as the onset point of the diamagnetic signal. Similar 
values were obtained from resistivity onset as well as ex- 
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trapolation technique used in Ref. [25| . 

The temperature scans of the dynamic susceptibiUty 
Rex at fixed dc fields from to 8T for the Y-1221 sin- 
gle crystal with H^^WHac along the c axis are shown in 
Fig.[T] Hereafter we deal only with Rex and omit "Re" 
for brevity; - this quantity is normalized so as to have 
— 1/477 in the limit Hdc and to correspond to the 
full expulsion of the ac signal for S small relative to the 
sample size. There are clear anomalies in high-field scans 
the locus of which appear as a continuous line. The inset 
to Fig.[T] shows magnetic field scans at constant T with 
a clearly visible anomaly. The positions {H, T) of the 
anomaly in field sweeps coincide with those observed in 
X {T) scans. The field Hc2{T) can also be extracted from 
the data as points where x first deviates from the upper 
envelop corresponding to the normal state skin effect; the 
observed anomaly clearly differs from Hc2- 
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FIG. 1: (Color online) x (T) in Y1221 measured in fixed fields 
with 1 kOe step. The anomaly corresponding to the VL tran- 
sition is clearly seen in the upper left corner (region of red 
color online). Inset: x (H) scans showing the same anomaly. 



Figure [5] shows x (T) for Y-1221 (upper panel) and 
Lu-1221 (lower panel) with a focus on low temperatures 
and high fields where the anomaly clearly seen in both 
samples. Hence, the observed behavior is not limited to 
a particular sample, but is seen on both mirror furnace 
grown Y-1221 and flux-grown Lu-1221. 

Figure [3] is a summary of the data in reduced coor- 
dinates T/Tc and iJ/i?c2 (0). The curves i?c2/i?c2(0) 
are nearly the same for two samples, although as shown 
in the inset, the absolute values are somewhat differ- 
ent. For comparison, we also show magnetization data 
of Ref. 27[ and the low-field branch of the R-S transi- 
tion according to SANS measurements of Ref. [l3| (shown 
with the uncertainty corresponding to various criteria for 
the peaks splitting that mark the R-S transition). The 
anomaly amplitude decreases with decreasing field and 
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FIG. 2: (Color online) Selected areas of x {T, H) curves show- 
ing VL transition in Y1221 (top panel) and Lul221 (bot- 
tom panel). At each fild the characteristic temperatures 
corrsponding to Hc2 (T) line as well as (T) are indicated 
by arrows at the lower frame. 



raising temperature (Fig. [J) , and the data cannot be ex- 
tended to the point where they are anticipated to merge 
with the low field and high temperature SANS results. 

The VL response to a small amplitude ac perturbation 
is determined by the length (23, 24| given by 



(1) 



where the term containing the Labusch pinning param- 
eter K is the Campbell's (in high fields, we do not 
distinguish between H and the magnetic induction). In 
the second formula Ac is expressed in terms of the disor- 
der length Vp and the critical current jc for the model of 
weak collective pinnin g ex pected to hold for high quality 
crystals of this work [23|. For weak disorder and small 
jc, A ^ Ac; also, one can use the relation — + (u^), 
where coherence length ^ is taken to be the range of the 
pinning potential and (u^) is the mean square thermal 
fluctuation of vortex positions . 

The VL shear moduli vanish in the limit H — > H^, the 
lattice softens, and (u^) diverges. As argued in the 
nonlocal corrections to the isotropic intervortex interac- 
tion responsible for the square structure are washed away, 
and the lattice near Hc2 is always triangular (for tetrago- 
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FIG. 3; (Color online) Phase diagram obtained from mea- 
surements of Hc2 (T) and Hq (T). Inset shows the same data 
in reduced units. Also shown are the results from Refs. [23] 
and [H. 



nal crystals in fields along c axis). With decreasing field, 
the fluctuations weaken, and at some field i?n < Hc2, the 
lattice acquires the square structure. The SANS data on 
cubic VsSi support this picture [18|. It has been esti- 
mated in Ref.ly that this happens when (u^) reduces to 
values on the order for clean materials. In most of 
the domain < H < Hc2, the fluctuations amplitude 
exceeds ^, and for a qualitative discussion we can take 



cH 



(2) 



The thermal average (u^) depends on all VL elastic 
moduli. The number of independent moduli and their 
values depend on the lattice structure. For flelds small 
on the Hc2 scale, the moduli were counted and evaluated 
in [2^ (and in [i^] within the Ginzburg-Landau scheme 
near T^). The point relevant for our discussion is that the 
number of independent moduli depends on the symme- 
try of the lattice (as within the standard elasticity theory 
for anisotropic solids [1^). When the VL rhombic unit 
cell with diagonals aloiig [100] and [010] of tetragonal 
borocarbides (see [1, [1, [13, [HI, [3, [11]) transforms to a 
square, it acquires an extra symmetry associated with 
the four-fold rotational axis. This reduces the number of 
independent shear moduli because the shears polarized 
along [100] and [010] (which are different in the rhom- 
bic lattice) become identical. In addition, the squash 
modulus for the deformation responsible for the transi- 
tion turns zero at the transition point [13, IS]- One, 



therefore, expects 
transition point. 



be different on two sides of the 



The calculation of this difference at the transition in 
high fields would have required evaluation of the elastic 
moduli in the presence of thermal fluctuations, a non- 
trivial problem to be addressed elsewhere. Still, we can 
state that, with increasing scattering, the domain of the 
square lattice should shrink and the R-S phase boundary 
should move to lower reduced fields H\j/Hc2 0- Physi- 
cally, this is because the transition to the square struc- 
ture is caused by the tetragonal symmetry of the Fermi 
surface and of the order parameter which enter the in- 
tervortex interaction through nonlocal corrections to the 
London theory [^ . In low fields, the corrections add the 
square symmetric interaction at distances on the order 
of the nonlocality range p which for clean superconduc- 
tors is on the order of foi zero temperature coherence 
length. This range, however, shrinks when the mean-free 
path becomes shorter. Hence, one has to go to higher 
fields for the square-symmetric interaction to cause the 
R-S transformation. In the high field domain with fluc- 
tuations induced lattice disorder, one has to go to lower 
fields with weaker fluctuations to reach the domain with 
sufficiently strong nonlocal effects to cause the R-S tran- 
sition. This happens when p ~ (w^). Since p is reduced 
by scattering, the transition takes place at smaller (it^), 
i.e., at lower fields. 
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FIG. 4: (Color online) Relative jump in x at the anomaly 
in YNi2B2C as a function of HHc2 (T). Solid line is a fit to 
Eq.(|9]). Inset: the definition of the anomaly size, AttAx- 



To estimate the dynamic susceptibility we note that in 
our data — 47rx ^ 0.1 at the location of the transition. 
This implies that in our case S is comparable or larger 
than the sample size. A weak screening of the ac field in 
a disc-shaped sample of a radius R is described by the 
London-like equation (see, e.g., [3]) A — (5^V^A — for 
the vector potential A = hxr/2-|-a, where a is a small 
correction to the uniform field. This gives 



h X r/2 : 



(3) 
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Given the current j, one readily obtains the magnetic 
moment and x = — -R^/327r5^. Hence, 



In the sample with 9% Co, 



Attx 



852 - 2cH ^ ' 



(4) 



If {v?) varies by /^{v?), the change in the measured 
quantity is: 



A(-4^X) = - 



(5) 



Since jc/{u^y^^ decreases on warming, one expects the 
anomaly size to decrease with raising T seen in our data. 
The relative change in the susceptibility is: 



Ax 

X 



A(u2) 



(6) 



At the transition (^2) ^ for clean materials. As is 
seen at Fig|4l Ax/x ^ 10^2 in large fields. Hence, our 
method is sensitive to small changes in thermal averages: 
yA(^~eo/10. 

The elastic properties of the incompressible square lat- 
tice in fixed fields are characterized by four independent 
moduli which can be chosen as the shear c^; along one of 
the square diagonals x = [100], the squash Cs transform- 
ing the square to a rhombus, the rotation Cr, and the tilt 
C44. At the transition Cs — 0, whereas according to 



Cx (1)0 Hc2{0) 



(7) 



(p varies from 0.67^o to 0.3^o on warming from to Tc 
in the clean limit). In this situation, Ref. 6] gives {u^) ~ 
f{H,T)ln{cx/cr), where / d oes not depend on elastic 
moduli. As is shown in [gI, l2d|, only the shear moduli cge 
change sharply at the transition, so that we estimate: 



A(m2 



AC66 



C66 In(c66/Cr) 



(8) 



The ratio Acge /cge should not depend strongly either on 
H or on T since the change of the shear moduli is of a 
geometric nature (due to the symmetry change) . We thus 
expect the following behavior of the measured quantity: 



Ax C 



(9) 



with C K, const. Fig. [4] shows the data Ax/x plot- 
ted against HHc2{T); the solid line is calculated tak- 
ing i?c2(0) = 8.1 T with the fit parameter C = 0.75. A 
qualitative nature of our argument notwithstanding, we 
consider this result as supporting our interpretation of 
the X anomaly as caused by the R-S structural trans- 
formation of the vortex lattice. It should be noted, 
that the qualitative model we use disregards pinning rel- 
ative to the thermal fluctuations: the disorder length 



this is unlikely to be true. This might be the reason why 
the anomaly Ax/x is nearly washed away in this sample. 
Large scale inhomogeneities may have the same effect. 



A similar claim has been made in Ref. 27| in which an 
unusual behavior of the reversible magnetization Mdc{T) 
was interpreted as caused by the R-S transition. The po- 
sition of the anomaly on the HT phase diagram reported 
here differs from that of Ref. [27| . We did not find an 
anomaly in Mdc{T, H) in our samples. 

To conclude, the tunnel-diode resonator is proven to 
have enough sensitivity to measure small changes in 
the rf penetration depth in the high-field mixed state 
(i.e. of the Campbell length) which accompany the VL 
structural transformations. The high-field branch of the 
rhombic-square lattice transition caused by thermal fiuc- 
tuation of vortices is mapped up. 
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